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We show that under a general disformal transformation the linear comoving curvature pertur¬ 
bation is not identically invariant, but is invariant on superhorizon scales for any theory that is 
disformally related to Horndeski’s theory. The difference between disformally related curvature per¬ 
turbations is found to be given in terms of the comoving density perturbation associated with a 
single canonical scalar field. In General Relativity it is well-known that this quantity vanishes on 
superhorizon scales through the Poisson equation that is obtained on combining the Hamiltonian 
and momentum constraints, and we confirm that a similar result holds for any theory that is dis¬ 
formally related to Horndeski’s scalar-tensor theory so long as the invertibility condition for the 
disformal transformation is satisfied. We also consider the curvature perturbation at full nonlinear 
order in the unitary gauge, and find that it is invariant under a general disformal transformation 
if we assume that an attractor regime has been reached. Finally, we also discuss the counting of 
degrees of freedom in theories disformally related to Horndeski’s. 


I. INTRODUCTION 

A primordial epoch of inflation and the late-time accelerated expansion of the universe constitute two key elements 
of the standard model of modern cosmology, which is in very good agreement with observational data. Many of the 
models proposed to try and explain these two epochs of accelerated expansion rely on the introduction of an additional 
scalar degree of freedom, either in the form of an unknown scalar field in the matter sector, such as an inflaton or 
quintessence field, or as part of a modified gravity sector, such as in f(R) gravity or Brans-Dicke scalar-tensor gravity. 

Recently, efforts have been made to determine the most general form of scalar-tensor theory that encompasses the 
examples mentioned above and more. In the spirit of effective field theories, such a theory would allow one to introduce 
a common parameterisation for a wide range of models, making it much easier to understand the relation between 
different models and to compare their predictions with observations. In trying to construct the most general form of 
scalar-tensor action, a key requirement is that the resulting equations of motion are second order in time derivatives. 
The appearance of higher order derivatives is generally associated with the presence of additional, ghost-like degrees 
of freedom, and the associated Hamiltonian becomes unbounded from below, both in the presence of even higher 
order derivative terms [1] and odd higher order derivative terms [2]. As a result, if the system is coupled to another 
“normal” system, then the total system will develop a so-called Ostrogradsky instability. The most general form of 
scalar-tensor action that gives rise to second order equations of motion was derived by Horndeski over 40 years ago 
[3], and was rederived just a few years ago in the context of so-called Galileon models [4-7]. More recently, however, 
it has become apparent that there exist theories that do not belong to Horndeski’s theory but that nevertheless do 
not suffer from Ostrogradsky instabilities, propagating only 3 degrees of freedom [8-15]. Whilst these theories may 
appear to give rise to higher-order equations of motion at the level of the Euler-Lagrange equations - which is why 
they are not included in Horndeski’s theory - it has been shown that the higher-order time derivatives can be removed 
by making use of the time derivative of a special linear combination of the gravitational equations of motion, thus 
rendering the equations of motion second order with respect to time derivatives [16]. These theories are therefore 
interesting and their phenomenology has been investigated in the literature, e.g. the screening mechanism [17-19] and 
possible observational signatures [20]. 

In exploring this class of general scalar-tensor theories, use is often made of disformal transformations of the metric, 
which take the form [21] 

9nv = «(</>, ^)0n v + /5(0, X)d li (j)d v (l), (1.1) 
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where X = —g fJ-u d /J ,(f>d v (j)/2. This is a generalisation of the more familiar conformal transformations, for which 
/3(4>,X) = 0 and a(4>,X ) —> a((f>). The different representations of a theory, written in terms of disformally related 
metrics, are often referred to as being written in different “frames”. In some cases a transformation of the form (1.1) 
can be used to remove non-minimal coupling between the scalar field <f> and the Ricci scalar or Einstein tensor at 
the level of the action, leaving only a canonical Einstein-Hilbert term [8, 22, 23]. This particular frame, if it exists, 
is referred to as the Einstein frame. It is known that the form of Horndeski’s action is preserved under disformal 
transformations if a and /? only depend on (j> [8, 23]. Similarly, the form of the theories beyond Horndeski’s considered 
in [9, 12] is known to be preserved under disformal transformations with a = a(cf>) and /3 = /3(cj),X). Allowing for 
an X-dependence of a, however, allows one to transform between theories belonging to the class considered in [9, 12] 
and those that lie outside it [8, 12]. 

In many cases it is easier to solve for the dynamics of the scalar field </> coupled to gravity if we first rewrite the 
theory in terms of a metric that is disformally related to the original metric as in (1.1). Having solved for </> and g 
however, it is often the case that we would like to relate these quantities back to the original metric g^ v , for example 
if matter is minimally coupled to this metric, i.e. it defines the so-called Jordan frame. In the context of cosmology 
we are particularly interested in the transformation properties of perturbations, and especially the so-called comoving 
curvature perturbation, 1Z C , which is defined as the curvature perturbation on time-slices of constant (f>. It has been 
known for some time that the comoving curvature perturbation is invariant under conformal transformations with 
/ 3(cf>,X) = 0 and a((f>,X) —> a(cf >), both at the linear level [24] and the fully nonlinear level [25, 26]. More recently, 
the invariance of 1Z C was also confirmed for transformations where a((j>,X) —> a(</>) and /?(</>, X) — > /3((j)) [27], and 
finally for the case where a((j>,X) —»■ and /3 has both <f >- and X-dependence [28]. The disformal invariance of 

the comoving curvature perturbation is a very useful result, as it means that if we are ultimately only interested in 
the comoving curvature perturbation then we are free to solve the system in whichever frame is most convenient. It 
is thus natural to ask whether or not the disformal invariance of 1Z C holds for the most general form of disformal 
transformation, where we allow for an X-dependence of both a and /3 , and this is the question we address in this paper. 
We will show that the case where an X-dependence of a is included is crucially different to the previously considered 
cases, and the comoving curvature perturbation is not identically invariant under such a disformal transformation. 

The paper is organised as follows. In §11, we investigate the transformation properties of linear perturbations, 
choosing to leave the gauge unfixed and to work with gauge-invariant quantities. We elucidate that the comoving cur¬ 
vature perturbation is not identically invariant under disformal transformations when one allows for an X-dependence 
of a , but that the difference is given in terms of the gauge-invariant comoving density perturbation associated with 
a single canonical scalar field. We then show that in the context of Horndeski’s theory a sufficient condition for 
the comoving density perturbation to vanish is 1Z C = 0, and that under reasonable assumptions - i.e. neglecting the 
so-called decaying mode of TZ C - this condition is satisfied on super-horizon scales. Details of this calculation are 
presented in Appendix A. Consequently, we conclude that on superhorizon scales, and under reasonable assumptions, 
the comoving curvature perturbation is disformally invariant for any theory that is disformally related to Horndeski’s 
theory. In §111, we consider the transformation properties of perturbations at the nonlinear level, and in this analysis 
we find it easier to make use of the unitary gauge, where 5<p = 0, which makes /3 irrelevant to the transformation 
law of the comoving curvature perturbation. We find that the comoving curvature perturbation is invariant under 
disformal transformations if we assume that an attractor regime has been reached. In such an attractor regime X can 
be re-expressed as a function of q i, meaning that the situation is exactly the same as in the case of an X-independent 
a. Using our result, we deduce that in the attractor regime the nonlinear curvature perturbation is conserved on 
superhorizon scales in any theory that is related to Horndeski’s by a general disformal transformation. In §IV, we 
discuss the counting of degrees of freedom in theories disformally related to Horndeski’s. Given that we know Horn¬ 
deski’s theory to be healthy, i.e. its equations of motion are second order and it propagates only 3 degrees of freedom, 
we focus our attention on how the counting of degrees of freedom is affected by a disformal transformation. Using 
a toy model we demonstrate that, as one would naively expect, the number of degrees of freedom as determined by 
a full Hamiltonian analysis should be unaffected by a disformal transformation, which is agreement with the recent 
analysis in [29]. §V is devoted to conclusions. 


II. LINEAR ANALYSIS 

We consider the disformal transformation given in (1.1). To make it a well-defined redefinition of fields, we require 
that the inverse disformal transformation from the tilded frame to un-tilded frame also exists, which requires the 
existences of g and the solvability of X = X(</>, X). 

Firstly, from (1.1) we can determine that the inverse metric g^ is given as 

g» v /?((£, X)^^ 


a^X) a(0,X)[a(0,X)-2X/3(0,X)] 


( 2 . 1 ) 
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where d^cj) = g^d^cj). From this expression we see that the inverse g will exist as long as 

a(a - 2X/3) ± 0. 


( 2 . 2 ) 


This could also have been inferred from the relation between the determinants of g^ and g^, which can be derived 
by contracting (1.1) with g^ and taking the determinant [30] 

g^~g„ a = a(s» a + ?-d^d a A => — = a 3 (a — 2Xf3), (2.3) 

V a ) g 


where g and g are the determinants of g and g^ respectively. 

Secondly, let us consider the solvability condition for X in terms of <fr and X. The complication that arises in trying 
to invert the transformation given in (1.1) is the appearance of g^ v in X. As discussed in [31], in order to be able to 
invert the transformation we thus need to express X = —g^d^d^cf )/2 in terms of X = —g^d^tfidvrf)/ 2. When (2.2) 
is satisfied, we can contract (2.1) with d^ifid „</> and obtain the relation 


X = 


X 

(*(</>, x)-2 xp&xy 


(2.4) 


The solubility of this relation for X requires dX/dX ^ 0, namely, 

a - lax + 2X 2 p x , 

(a - 2A/?) 2 ^ ’ 


(2.5) 


where ax = dxa and fix = dxP* 1 The same condition can also be determined by considering the Jacobian 
dg^v/dgap and requiring its determinant to be non-vanishing, as in [8]. Combined with (2.2), 


a(a - 2 Xp){a - Xa x + 2X 2 f3 x ) ± 0 


( 2 . 6 ) 


is the necessary and sufficient condition for the invertibility of the disformal transformation. If (2.6) is satisfied, we 
can solve (2.4) for X = X(<f>,X) and obtain the inverse disformal transformation 


1 


9n* = 


~9iif 


Pi&X^X)) 


d^dvCp, 


&(<!>, X(<j>,X)) a(</>,X((f>,X)) 

In the following, we consider disformal transformations that satisfy the invertibility condition (2.6). 
We first focus on a linear analysis of perturbations. Let us take a line element of the form: 

ds 2 = g^dx^dx 1 ' 

= — (1 + 2 A)dt 2 + 2a(diB — Si)dx l dt 

+ a 2 [(l — 2 + 2 didjE + diFj + djFi + hij]dx l dx ^, 


(2.7) 


( 2 . 8 ) 


where F) and satisfy F iti = ha = hijp = 0.* 2 
ds 2 = g^ u dx ti dx v 


We then consider a disformal transformation of the form (1.1). 


— [— 00(1 + 2A) — 8a + (/3o + 5/3)<(> 2 + 2/3o<( , <5<)6]dt 2 
+ 2[aoa{diB — Si) + doi’diS^jdtdx 1 


CM.QCL 


r \ 

1 — 2y/i -|—— ) Sij + 2 didjE + diFj + djFi + hij 

o-o J 


dx l dx\ 


(2.9) 


where a dot denotes d/dt and we have decomposed a = a$+ 5a and (3 = do + 8/3. Note that we are choosing to leave 
the gauge unfixed, which will make the interpretation later on more obvious. We would then like to rewrite this new 
line element in the same form as (2.8), namely as 

ds 2 = — (1 + 2 A)dP + 2 a{diB — Si)dx l dt 

+ a 2 [(l — 2 %j)5ij + 2didjE + diFj + djFi + h i j]dx l dx : ’. (2-10) 


" 1 We will adopt a similar notation for derivatives with respect to <b. e.g. ap = dpa. 
2 Note that repeated indices are summed over. 
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At background level this then gives us 


= a\fa 0 and dt = \J a o — Pop 2 dt 


( 2 . 11 ) 


where we note that a dot still corresponds to taking the derivative with respect to t rather than t. At the level of 
perturbations, from the 00-component we obtain 


A = 


- l —r- (a 0 A + ]-5a - \<j> 2 5(] - p 0 p5p ) . 

a 0 - Pop \ 


From the 0*-component we have 

B = 


- B + — 50 

a 0 - pop 2 V a ° a ) 


and Si = Si. 


Finally, from the *j-component we find 


0 = 0 — 


5a 

W 


( 2 . 12 ) 


(2.13) 


(2.14) 


whilst E, Fi , and hij remain unchanged. As such, we see that the vector and tensor perturbation are invariant under 
a general disformal transformation at linear level. Whilst the above expressions are problematic if ao — Pop 2 vanishes, 
recall that we are assuming a — 2X/3 ^ 0 for the existence of g 111 ', which yields «o ~ Pop 2 ^ 0 at background level. 3 
Turning to the gauge-invariant comoving curvature perturbation, 1Z C , it is defined in the original frame as 


Kc = -0 - -50, 

9 

where H = a/a. In the new frame we similarly have 


Tl c = —0 — 


H 




d(f>/dt 

where H = (1 /a)da/dt. Using the background relations (2.11) we have 

1 


d<f> 

dt 


and H = 


ao — 


on ^ PoP 2 


H + 


ao 

2ao 


(2.15) 


(2.16) 


(2.17) 


which, on combining with (2.14), gives us 


K c = Tic + (5a - ^ 

2 ao 


(2.18) 


In the case that a = a(<j>) we have 5a — do 5p/p = 0, meaning that 7 l c = 7 l c , which is consistent with [27, 28]. 
However, in the case that we allow for an A-dependence of a , we more generally get 


5a —3(50 = aox \ 5X — 


(2.19) 


where a 0 x = dao/dX and X 0 = 0 2 /2. Using 5X = <fi(5p — pA), we thus find 

A -a — a ° x 
Tl c Tic — „ £s ; 

2a 0 


( 2 . 20 ) 


* 3 In fact, the requirement that the perturbative expansion be valid in the disformally related frame, e.g. A <C 1, will put tighter constraints 
on the background-dependent coefficients appearing in eqs. (2.12)-(2.14) than those imposed by the invertibility of the transformation, 
but we will not consider this issue here. 
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where 


e s = <p{Sp — pA ) — pSp = SX — pSp. 


( 2 . 21 ) 


We thus see that the comoving curvature perturbation is not identically invariant under disformal transformations 
with a = a{p : X). Note that if one takes the gauge Sp = 0, 1Z C as determined by (2.20) coincides with the quantity 
(new defined in Eq. (96) of [15]. * 4 The importance of Cnew was discussed in [15]: it absorbs all terms in the action 
generated by a disformal transformation that explicitly depend on the time derivative of the perturbation of the lapse 
function. In so doing, it makes it explicitly clear that no additional degrees of freedom appear as a result of the 
disformal transformation. In light of the above analysis, we see that the appearance of the quantity £new(= 1Z C ) is in 
fact very natural, as it simply corresponds to the transformed comoving curvature perturbation. 

The quantity e s is a gauge-invariant quantity corresponding to the gauge-invariant perturbation of X, as seen in 
(2.21). It also coincides with the comoving density perturbation for a single canonical scalar field, e s = Sp s = Sp—3HSq, 
where Sp = SX + V^Sp is the density perturbation and 8 q = —pSp is the velocity potential for the energy momentum 
tensor of the scalar field [32] . e s is also related to the intrinsic entropy perturbation of a canonical scalar field as 


f Sp _ < 5 p\ __214_ £ 

\J~JJ~ 3<p 2 (3Hp + 2V^) €s ' 


( 2 . 22 ) 


So far, we have not assumed any particular scalar-tensor theory, and thus (2.20) holds for any theory. We now 
proceed to consider specific theories, in order to determine how e s behaves. In General Relativity, if the scalar field 
is the dominant energy component of the universe, then from Einstein’s equations we are able to determine that e s 
satisfies the Poisson equation 


--T = ' 


2 ’ 


(2.23) 


where ’P = ip + a 2 H(E— B/a) is the gauge-invariant Bardeen potential. As such, e s is suppressed by k 2 on large scales 
as long as 'k remains finite, which implies that the difference between 1Z C and 1Z C will also vanish on large scales. 

It has also been shown that in a subclass of Horndeski’s scalar-tensor theory e s still vanishes on superhorizon scales 
if the comoving curvature perturbation remains constant [27]. To the best of our knowledge, however, it has not yet 
been explicitly shown for the full Horndeski theory, and this is what we will now proceed to confirm. Here we simply 


give the result, and more details can be found in Appendix A. 

Horndeski’s action takes the form 

S = j d 4 Xy/^C H , (2.24) 

with Ch = X)i =2 & an d 

C 2 = K(P,X), (2.25) 

C 3 = -G 3 (P,X)DP, (2.26) 

U = G 4 (0, X)R + Gax [{Op) 2 - (V M V^)(V^V»] , (2.27) 


£5 = G 5 ( 0 , X)G^WX V P - i G 5X [(D 0) 3 - 3 {np)(X^„p){V»V v p) + 2{X^X a p)(X a V^XV^V^)] . (2.28) 

Focusing on scalar perturbations, and taking the spatial gauge E = 0 in (2.8), the equations of motion for A, B , p 
and Sp were derived in [33]. Combining the constraint equations that result from varying the second order action 
with respect to A and B, we are able to derive the following Poisson equation for e s 


e s 


k 2 ftCxH 2 
a 2 H 2 CxA 4 - C 3 A 1 


AilZ r 


— (T 
H K 


TIP) 


(2.29) 


where the coefficients A, and Ci depend only on background quantities and are given in Appendix A. In the case of 
General Relativity with a single canonical scalar field, where G 3 = G 5 = 0, K{p,X) = X — V(p) and G 4 = 1/2, we 
have 


A 1 = 6 tf, -A 3 = 2, A 4 = 2X- 6H 2 , A 5 = -2H, C x = 2 and C 3 = -2 H, (2.30) 


* 4 In fact, the two expressions do not exactly coincide, but this is due to a typo in Eq. (96) of [15], where N should be replaced by 
N(= 1 + SN), so that to linear order the expression for £ n ew should be £new = £ H— ^-SN. We thank J. Gleyzes for confirming this 
point. 





6 


so that (2.29) reduces to (2.23). 

Given the form of (2.29), one can conclude that e s vanishes on superhorizon scales i.e. in the limit k <C aH 
so long as the coefficient of k 2 /(aH ) 2 on the right hand side is finite in this limit. Alternatively, from the momentum 
constraint we have 

e s = (2.31) 

G 3 

which is given in (All) but we repeat here for convenience. As such, we see that a sufficient condition for the vanishing 
of e s - and thus disformal invariance of 1Z C - is that 1Z C is conserved. 

In the case of General Relativity plus canonical scalar field, it is well known that 1Z C is conserved on superhorizon 
scales, provided that the so-called decaying mode can be neglected. Explicitly, one finds that 1Z C oc l/(ea 3 ), where 
e = (j) 2 /{2H 2 ). As such, we see that 1Z C is indeed decaying — and therefore negligible — provided e is not decaying 
faster than a~ 3 . This condition is satisfied in almost all standard slow-roll inflation models, where the slow-roll 
parameter e is itself taken to be slowly varying. There is, however, a special class of inflation models - dubbed 
“ultra-slow-roll” models - for which extra care is needed [34-37]. In the simplest ultra-slow-roll inflation model with 
constant potential one finds e oc a -6 , which means that 7 Z c is growing as a 3 . Interestingly, however, we still find that 
e s decays as a -3 , which follows from the fact that e s and 1Z C are related by a factor of </> 2 oc a~ 6 . As such, provided 
the disformal transformation is such that the factor a^x/^o appearing in (2.20) is not growing faster than a 3 , we see 
that 7 Z c is disformally invariant even in the case of ultra-slow-roll inflation. 

In the more general case of Horndeski’s theory we have a similar result. As was shown in [7], on superhorizon scales 
TZ C has the two independent solutions 5 


/ 1 

1Z C = const. and 1Z C oc / —— ndt', 

J Qsa 3 


(2.32) 


where Gs = @1/Q 2 )Gt + 3Gt an d £> © and Gt are as defined in Appendix A. We thus find that 1Z C oc 1 /(Gsa 3 ), which 
is decaying provided Qs is not decaying faster than a~ 3 . In the case of standard slow-roll inflation we expect this to 
be the case, but there will be exceptions analogous to ultra-slow-roll inflation. Strictly speaking, even if the decaying 
mode can be neglected, in order to then conclude that 1Z C = 7 Z c on superhorizon scales we additionally must assume 
that the combination aoxC\(j > 2 / { 010 C 3 ) is not growing faster than 1Z C is decaying. This again seems reasonable if we 
assume that background quantities are evolving slowly, but perhaps there may be some exceptions in the very general 
context of Horndeski’s theory. Conversely, as was the case with ultra-slow-roll inflation, even if 1Z C is not decaying, 
1Z C may still be disformally invariant if the combination aoxC\(jr /is decaying faster than 1Z C is growing. 

To reiterate, our main conclusion is that in the context of Horndeski’s theory, a sufficient condition for the vanishing 
of e s - and thus disformal invariance of 1Z C - is that 1Z C is conserved, and this is the case on superhorizon scales so 
long as we can neglect the so-called decaying mode of 1Z C . Models in which the decaying mode cannot be neglected 
- such as ultra-slow-roll inflation - must be considered on a case-by-case basis, but interestingly it seems that the 
disformal invariance of 1Z C does not necessarily break down in such cases. For the remainder of this section we will 
restrict ourselves to considering models in which the decaying mode can be neglected. 

Using the above results, we can argue that the comoving density perturbation e s should vanish on superhorizon 
scales in any theory that is disformally related to Horndeski’s theory as follows. Suppose we have two theories, theory 
A and theory B, that are both disformally related to an element of Horndesdki’s theory, theory H. The metrics for 
these theories are related as 


9 $ = a (H ^ A) 9^ } + 

9 ^=a^ 9 ^ + ^d^. 

We can then consider a disformal transformation between theory A and theory B 

9= Oi {A ^ B) g (A) + /? (A ^ B) d^d v (t>, 

with 


v (a^b) _ 


y(H-)-B) 


q,(H—>A) > 
^j(A—>B) = ^j(H—>B) _ 


™(H^B) 

(A__fl(H^A) 

«( H ->A) P 


(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 


*5 


Note that these two solutions are in fact valid on scales larger than the sound horizon of the scalar perturbation 1Z C . Horizon crossing 
is defined by c^k 2 = a 2 H 2 , where c s is the sound speed of 7 Z c and is in general different from unity. See [7] for the general expression. 
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Then, the comoving curvature perturbations in these theories are related as 


- n ( H) 


7^ b) - 7^ h) 


(H^ A) 

* 0 A 


f (H) 


2 a, 


(H-vA) s 


(H-s-B) 

x 0X 


,(H) 


2 a, 


o 


(2.38) 

(2.39) 


Here, is understood as a derivative with respect to = g'^d^cpdutp. As the comoving density perturbation 

ei H ' ) vanishes on superhorizon scales in Horndeski’s theory, lzi A ^ = IZ^ = 7Z^ on superhorizon scales. By considering 
the disformal transformation of 7 Z c between theories A and B, which gives us 


- izi A) 


(A—>B) 

Q 0A' (A) 

2 a^ B) 8 ’ 


(2.40) 


(A) 

the vanishing of the left hand side on superhorizon scales allows us to infer the vanishing of e« ' on superhorizon 
scales. 

Indeed, we can also confirm the above statement explicitly by considering how e s transforms under a disformal 
transformation. Using (2.12), (2.17) and (2.21), we obtain 


e s 


ao — «ox 2 + 


2/3oa (f)' 


(a 0 - M 2 ) 2 


(2.41) 


which recovers the result in [27] when a and /3 are functions of <f> only. Interestingly, even for a general disformal 
transformation with X-dependent a and /?, e s is disformally invariant up to a coefficient depending on background 
quantities. This implies that any theory disformally related to Horndeski’s theory should also have vanishing e s on 
large scales. 

Let us recall here that we are assuming the invertibility condition (2.6) is satisfied. This condition precisely guaran¬ 
tees that the coefficient in front of e s in (2.41) neither diverges nor vanishes. As such, the conclusion that e s vanishes 
if e s vanishes holds for any disformal transformation that is invertible. It is interesting to note that the coefficient on 
the right-hand side of (2.41) exactly coincides with the left-hand side of (2.5) evaluated at background level, i.e. it 
coincides with dX/dX. This can be expected, however, as e s corresponds to the gauge-invariant perturbation of X 
and e s to the gauge-invariant perturbation of X (recall (2.21)). 


III. NONLINEAR ANALYSIS 

In going beyond linear perturbations let us take the unitary gauge from the outset, where 8 (j) = 0. In this case, the 
spatial part of the metric is not affected by the disformal part of the disformal transformation, /3, due to the fact that 
di<f> = 0. Let us take the metric with nonlinear perturbations as 

ds 2 = -N 2 dt 2 + a 2 e 2 Kc 7ij ( dx i + A r*dt) {dx j + Wdt) , (3.1) 


where 


7^- =e 2did j E+d i F :j +d j F i +h i: j ^g2) 

and Fi and hij once again satisfy F l: i = ha = hij t i = 0. Note also that N l contains both scalar and vector components, 
and that spatial indices should be raised and lowered with y 1 - 3 and 7 ij, respectively. If we further parameterise a(4>, X) 
as a = aoe 2A “, where ao corresponds to the background part of a and e 2A “ contains nonlinear perturbations from 
this background value, then we see that in the gauge <5<^> = 0 the metric transforms as 


g 00 = -aN 2 + pft + aoa 2 e 2n c +2Aa 7 io N l N\ 
g 0i = aoa 2 e 2 ' R ' c+2Aa 'YijN 3 , 


Qij — aoct, e 


2 0 27J C +2A a 


7 ij, 


(3.3) 






from which we deduce 


N 2 = aN 2 — /3(f) 2 , a 2 = aga 2 , TZ C = 1Z C + Act, 


(3.4) 


whilst E , Fi , N l and hij remain unchanged at nonlinear level, which is consistent with the result at linear level in the 
previous section. As such, the vector and tensor perturbations are invariant at the nonlinear level. 

In analysing the last relation in (3.4), note that as /3 is irrelevant for the transformation law of 1Z C , the situation is 
equivalent to determining how 1Z C transforms under a conformal transformation. As such, in the case that a = a(<j)) 
our conclusion is the same as that reached in [26]: in the unitary gauge 5(j) = 0 so that Aa = 0, meaning that 1Z C is 
invariant. This is also in agreement with [28] and - at the linear level with the results of §11 and [27]. 

In contrast, when we allow for an A-dependence of a , even in the unitary gauge we find that Aa ^ 0 as a result 
of the dependence of A on N. Explicitly, in the unitary gauge we have a = a(</>, cp 2 /(2N 2 )). As such, we see that a 
will only coincide with its background value, hence giving Aa = 0 and 1Z C = 1Z C1 when the perturbation of the lapse 
function vanishes. At linear order, this condition corresponds to A = 0, which is thus consistent with the requirement 
found in §11 that e s must vanish if we are to have TZ C = 1Z C , as in the unitary gauge we have e s = — <j) 2 A . 

In order to aid an intuitive understanding of this condition, let us define the proper time t as dr = Ndt. Starting 
with the definition of e s at linear order, we can see that it can be rewritten as 


e s 


d T 6d 2 T 0 


(S(d T </>) 
1 d 2 A 



(3.5) 


where the term in brackets corresponds to the relative entropy perturbation between d T <f> and </>. We can thus see 
that e s vanishes - in turn giving 1Z C = 7 Z c - when d T (j) = /(</>), where is some function of </>. Similarly, turning 
to the nonlinear case, we see that in terms of r we can write a = a(</>, (d T (f>) 2 / 2). Imposing d T (f> = f((f>) means that 
in the unitary gauge a is equal to its background value, which in turn gives us Aa = 0 and 1Z C = 1Z C . The condition 
d T (j) = /(</>) is familiar to us as the condition for an attractor regime (see e.g. [38]), and we thus conclude that in the 
unitary gauge and an attractor regime the curvature perturbation is disformally invariant at the nonlinear level. Note 
that the requirement to be in an attractor regime is not as restrictive as it may sound. Indeed, the vast majority of 
standard inflationary models satisfy this condition. 

It has been shown in [39] that even at the nonlinear level the comoving curvature perturbation has a mode that 
remains constant on superhorizon scales in Horndeski’s theory. Provided that the other so-called decaying mode can 
be neglected, this allows us to conclude that on superhorizon scales the nonlinear comoving curvature perturbation 
is both conserved and disformally invariant in any theory that is disformally related to Horndeski’s if one is in the 
attractor regime. 


IV. THE NUMBER OF DEGREES OF FREEDOM IN THEORIES DISFORMALLY RELATED TO 

HORNDESKI’S 


In the preceding sections we have investigated the transformation properties of the comoving curvature perturbation 
both at the linear and nonlinear level, with our conclusions holding for any scalar-tensor theory that is disformally 
related to Horndeski’s theory. In this section we discuss whether or not such theories are well behaved in the sense that 
they do not suffer from so-called Ostrogradsky instabilities. Such instabilities generically arise in theories possessing 
equations of motion that are higher than second order in time derivatives, indicating the presence of additional ghost¬ 
like degrees of freedom. As we will see below, the equations of motion derived from theories disformally related to 
Horndeski’s theory do seemingly contain higher-order derivatives, thus suggesting that they are unhealthy. However, 
this is somewhat at odds with our expectation given that the theories are simply related to instability-free Horndeski’s 
theory by a redefinition of fields. In the following we try to address this apparent contradiction. 

The action associated with the aforementioned class of theories can be written in the following form 

S = J d^x^J~^C H {g^A) + J d^Xy/^gCmtg^), (4.1) 

where we have Horndeski’s Lagrangian, Ch , written in terms of the metric g^ and matter is minimally coupled to 
the metric g^. The metric is disformally related to g IJV as in (1.1). The frame defined by g^ is the “Horndeski 
frame”, in which the gravitational Lagrangian coincides with that of Horndeski’s theory and matter is non-minimally 
coupled to the scalar field through g^ = ( g^ — pd^cjid^cj))/a. However, we are interested in the equations of motion 
in the Jordan frame, defined by g llVl in which matter is not coupled to the scalar field. Varying the above action with 
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respect to g pl , and 0 yields their equations of motion, 


a£ H + 2 £ h (axg P o 


,8xd p 0d a (j>) 8^08''0 - 


2uyja{a - 2X(3) 


rr = 0 , 


asja(a - 2Xp) {£% {a x g p<J + P x d p 0d <7 0) 8»0 - 2p£%'8„0} 


+ a^a(a - 2X/3) E^ (a^gpa + P < j,d p 0d rT 0) + E^ = 0, 


(4.2) 


(4.3) 


where and £^ are determined by varying Horndeski’s Lagrangian with respect to g pl/ and 0 [7, 39], namely, 
HV~y£n) = \J~-g{£ l H Sg pu + £^60), and = ~^= denotes the energy-momentum tensor associated 

with the matter Lagrangian. Note that here 8^0 = g^datj). As £^ and E^ are the equations of motion for 
Horndeski’s theory, they contain at most second order derivatives of 0 and < 7 ^, meaning that they also contain at 
most second order derivatives of g pi/ . However, as g pv contains X = — g^ v 8p,08 v 0/2, in principle it is possible for 
£'x and £^ to contain up to third order derivatives of 0. As such, we see that (4.2) will contain up to second order 
derivatives of but potentially third order derivatives of 0. Similarly, (4.3) will contain up to third order derivatives 
of g pl , and fourth order derivatives of 0 . 

As already mentioned, it is reasonable to expect that the appearance of seemingly dangerous higher order derivatives 
in (4.2) and (4.3) may be spurious, as we know that our theory is related to a healthy theory by a field redefinition. 
Indeed, following the analysis of [31, 40], it is possible to show that the gravitational equations of motion (4.2) are 
equivalent to the equations of motion obtained by varying the action with respect to g pv so long as the transformation 
(1.1) is invertible. Explicitly, we obtain (see Appendix B for details) 


CM v 
°H 



(4.4) 


where T anc j j s gj ven explicitly in terms of as 


rpIlU _ 


a 2 ^a(a-2XP) 


rjnfll/ _ 


d^ 0 d ' / 0 


2 (a + 2X 2 p x - Xax) 


T m (axgpa + Pxdp0d a 0) 


(4.5) 


Note that contains at most first order derivatives of 0. On substituting the above expression for £^ into (4.3), 
we obtain 


avMa - 2XP) 4 0) + 


r TP a 

m 

2 


2 X 2 a x p(p — ct<i>[ot + 2X 2 p x ) 
a (a + 2X 2 p x - Xa x ) 


—T^d 0 _-_ 2 ~ Y/3 _ 

a m 2a(a + 2X 2 p x -Xa x ) 

2 


d p 0d a 0 

' 9pa a + 2X 2 p x 


(a x g P * + Pxdp0d a 0)T£d»0 

Xfix[2X 2 axfi(t> ~ + 2X 2 p x )] \ 

a(a + 2 X 2 P-Xa x ) ) 


(4.6) 


= 0. 


If T0 r 0 does not contain any second order derivatives of g pu , then we see that (4.4) and (4.6) both contain at most 
second order derivatives of g pu . Moreover, they both contain at most third order derivatives of 0, which result from 

the second order derivatives of g pL , appearing in E^ and E^f 1 . Note that, as a result of the kinetic mixing between 
the scalar field and the metric in (4.3), the equation of motion (4.6) contains mixing terms between 0 and the matter 
energy-momentum tensor and its derivatives. The rich structure of kinetic mixing found in general scalar-tensor 
theories has many interesting consequences, as discussed in [41]. 

Although the situation has been marginally improved in moving from (4.2) and (4.3) to (4.4) and (4.6), we still have 
third order derivatives of 0 appearing, which would normally indicate the presence of dangerous additional degrees of 
freedom. One can superficially remove these third order derivatives by using the trace of the gravitational equations 
(4.4) - sometimes referred to as a “hidden constraint” - to find an expression for 0 in terms of lower order derivatives 
of 0 and g pill which can then be substituted into (4.4) [ 8 , 12]. However, as recently highlighted in [16], this does not 
offer a rigorous proof that the number of degrees of freedom is unaffected by the disformal transformation, as the 
trace equation itself is a dynamical equation that must still be satisfied. 

A more rigorous approach, which is in the same vein as that mentioned above and has been demonstrated in [16], 
is to find a suitable linear combination of the gravitational equations of motion (4.4) that does not contain 0 and 
allows you to express 0 in terms of at most first order time derivatives of 0 and g pv . On taking the derivative of 
this relation we are then able to find an expression for 0 in terms of at most second order time derivatives of 0 and 
gp.v, which can be substituted back into the original equations. The difference here is that the special combination 
of gravitational equations used does not itself involve third order derivatives of 0. In [16], Deffayet et al. showed 
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that such a combination does indeed exist for a very general class of models in the context of covariantized Galilcons. 
More specifically, they considered theories that in four spacetime dimensions are equivalent to (2.24) but with the 
replacements Gix{(t>,X) F^{(j>,X) and G^x{(^,X) —> F$((f>,X), where F 4 and F 5 are arbitrary functions of <j> and 
X. When evaluated in the unitary gauge, this Lagrangian is equivalent to the Lagrangian of the theories beyond 
Horndeski’s considered in [9], and as mentioned in the introduction, the structure of this Lagrangian is preserved 
under disformal transformations with a = and = /3(<f>,X). The models under consideration here, however, 
are related to Horndeski’s theory by a more general disformal transformation with a = a(<p,X) and /3 = /3(cj),X), 
meaning that they do not belong to the class analysed by Deffayet et al. As such, the appropriate combination of 
gravitational equations with which we can express <j) in terms of first order time derivatives of and g remains to 
be found. 

Another rigorous approach to determine the number of degrees of freedom of the theory is to perform a Hamiltonian 
analysis. Such an analysis has been performed for the newly-discovered theories beyond Horndeski’s in Refs. [11, 12, 
14], but in these analyses the unitary gauge was taken from the outset, which could affect the outcome of the degrees- 
of-freedom counting [16]. A Hamiltonian analysis that does not rely on fixing the gauge has been performed for an 
example Lagrangian in [16], but its generalisation is yet to be carried out. 

In the context of theories that are disformally related to Horndeski’s, we are particularly interested in how the 
Hamiltonian analysis of a theory is affected by a disformal transformation. In the case of a conformal transformation, 
where /3 = 0 and a = it has been demonstrated in the context of f(R) gravity that the transformation simply 
gives rise to a canonical transformation of variables, thus rendering the Hamiltonian analysis unchanged [42]. More 
general metric transformations - which include disformal transformations of the form ( 1 . 1 ) - were also considered 
recently in [29], where they reached a similar conclusion; we will comment further on their results shortly. Whilst a 
complete analysis of how the Hamiltonian analysis of a theory is affected by disformal transformations is beyond the 
scope of this paper, it is nevertheless possible to see why we might expect the degrees-of-freedom counting to remain 
unchanged. 

The key distinguishing feature of general disformal transformations taking the form (1.1) as opposed to conformal 
transformations with /3 = 0 and a = a((f>) is that they induce higher order derivatives of some of the fields in the 
theory at the level of the action. 6 However, the fact that these higher order derivatives always appear in a certain 
combination ■ as we will see with the help of a toy model below - means that they are always associated with additional 
primary constraints in the Hamiltonian analysis. If these constraints are first class, or second class and give rise to 
secondary second class constraints, then they will remove the spurious degrees of freedom associated with the higher 
derivatives appearing in the action as a result of the transformation. Indeed, this is the case in the example model 
considered in Sec. Ill of [29], where a derivative of the lapse, N , induced by the disformal transformation is always 
associated with a derivative of the spatial metric, 7 y, leading to a primary constraint involving the corresponding 
canonical momenta ttn and 7 r 1 - 7 . The “exorcising of Ostrogradsky’s ghost” with constraints was also discussed in [43] . 

In order to demonstrate the above idea, let us consider a simple toy model consisting of two degrees of freedom 
x(t) and z(t), with Cq = Cq(x,x, z, z). Here we have in mind that x represents the metric degrees of freedom g^ u 
while z the scalar field (f>. Assuming this Lagrangian to be regular means that we have two degrees of freedom and 
require four initial conditions. Analogous to a disformal transformation, we then consider that x = x(y,z,z), where 
y(t) represents the disformally related metric and the transformation depends on the analogues of <f> and X , z 
and z respectively. In terms of this new set of variables we have Cq = Co{y, y , z, z, z), i.e. due to the appearance of z 
in the transformation we have picked up a dependence of Cq on z. Following the constrained Ostrogradsky approach 
to systems with higher order derivatives - see e.g. [44, 45] - we introduce an auxiliary degree of freedom as w = z, 
which we impose by adding a Lagrange multiplier to the Lagrangian, i.e. we have 

C = C 0 (y,y,z,w,w) + X(w - z). (4.7) 


This Lagrangian depends on four degrees of freedom and contains up to first order derivatives, so in general eight 
initial conditions are required. On calculating the canonical momenta we obtain 


Py = 


dCo 

dy ’ 


Pz = -A, 


Pw = 


9£o 

dw 


and 


P\ = 0, 


(4.8) 


meaning that we have two obvious primary constraints, namely 


$1 = p z + A 0 and $2 = Px ~ 0. 


(4.9) 


To be more specific, we are interested in the case where both a and j3 have an X-dependence. As discussed previously, in the case where 
a = ot(4>) and (3 = (3((f)) Horndeski’s theory is mapped onto itself, and we already know that Horndeski’s theory propagates only three 
degrees of freedom. Similarly, in the case where a = a((f>) and (3 = (3((f), X), Horndeski’s theory will map onto a subclass of the theories 
introduced in [9], which were shown to be healthy without making use of the unitary gauge in [16]. 
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Note that these two constraints result from the introduction of the auxiliary degree of freedom, and in this sense 
do not correspond to the “additional primary constraints” referred to above. In order to see how these additional 
primary constraints appear we note that 

dCo dCo dx dCo dx 8 Cq OCq 3Cq dx 3 Cq dx 

Py dy dx dy dx dy &n Pw dw d'z dx dz dx dz' 


so that we find the additional primary constraint 

$3 =P V - F(y, z, w)p w « 0, 


where 


F(y, z, w) 


dx 1 
dy dx/dz 


(4.11) 


On constructing the total Hamiltonian and requiring that the constraints be preserved under time evolution we 
find an additional secondary constraint, which we label $ 4 . Requiring $4 to be preserved under time evolution 
does not give rise to any additional secondary constraints, and we find that all four constraints are second class. 
As such, although we started with four degrees of freedom requiring eight initial conditions, we found four second 
class constraints, meaning that we only have to specify four initial conditions, corresponding to only two degrees of 
freedom. This is in agreement with the fact that we had two degrees of freedom present in the original Lagrangian 
Co(x,x,z,z). 

The key point in the above analysis was the appearance of the third primary constraint $ 3 . As alluded to earlier, 
this arose due to the fact that z only appears in the Lagrangian in a specific combination with y , which is a consequence 
of the form of the analogue of the disformal transformation x = x(y , z, z). 

In fact, the toy model considered above falls into the more general class of models and transformations considered 
in [29]. To see this we introduce two auxiliary degrees of freedom, which allows us to re-write the Lagrangian in a 
form that is linear in velocities, i.e. we write 


C = Cq{x, u, z, w) + Ai (u — x) + A 2 (w — z). (4.12) 

When written in terms of w, the transformation no longer depends on the velocity of any of the fields, i.e. we have 
x = x(y,z,w). This form of transformation and the form of action given in (4.12) are then of the form considered 
by Domenech et al. in [29]. Assuming that the transformation x = x(y,z,w) is invertible, they show that such a 
transformation is a canonical transformation, so that the set of constraints and constraint algebra are left unchanged, 
which in turn means that the number of degrees of freedom remains the same. Note that the analysis of Domenech 
et al. does not rely on the original Lagrangian being regular, and so offers a general proof that the number of degrees 
of freedom is unaffected on making a disformal transformation, so long as the transformation is invertible. We thus 
conclude that theories disformally related to Horndeski’s theory via an invertible disformal transformation of the form 
( 1 . 1 ) are healthy, in the sense that they propagate only three degrees of freedom. 


V. CONCLUSION 

In this paper we have examined how the comoving curvature perturbation transforms under the general disformal 
transformations of the metric given in (1.1). We began by considering linear perturbations, and showed that whilst 
the vector and tensor perturbation are invariant, the comoving curvature perturbation is not identically invariant 
under a general disformal transformation. The difference between disformally related curvature perturbations is given 
in ( 2 . 20 ), and is written in terms of the gauge-invariant comoving density perturbation e s associated with a single 
canonical scalar field. In the context of Horndeski’s theory we used the Hamiltonian and momentum constraints to 
derive a Poisson equation for e s , and found that a sufficient condition for the vanishing of e s is that 1Z C is conserved, 
which is the case on superhorizon scales so long as the so-called decaying mode of 1Z C can be neglected. As such, we 
concluded that the comoving curvature perturbation is disformally invariant on superhorizon scales for any theory 
that is disformally related to Horndeski’s theory, provided that the decaying mode of 1Z C can be neglected. Using this 
result, we saw that under the same mild assumption the comoving density perturbation is suppressed on superhorizon 
scales for any theory that is disformally related to Horndeski’s theory, which we were also able to derive explicitly 
from the transformation rule for e s under disformal transformations. The relation between e s and e s , given in (2.41), 
shows that e s is equal to e s up to a proportionality coefficient that depends on background quantities. The coefficient 
is finite so long as the condition given in (2.6) for the invertibility of the disformal transformation is satisfied. Based 
on these findings, we conclude that in most cases of interest in the context of inflation, we are free to work in any 
disformally related frame when we wish to calculate the superhorizon curvature and tensor perturbations that are 
required in making predictions for inflationary observables. This is an extension of the conformal invariance of the 
curvature and tensor perturbations that is frequently exploited in the context of simple scalar-tensor theories, where it 
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is often much easier to perform calculations in one frame than in another. One can therefore expect that the disformal 
invariance of perturbations will also prove to be very useful in this wider class of theories. 

Using the unitary gauge, we also considered the comoving curvature perturbation at full nonlinear order, and found 
that it is invariant under disformal transformations if we assume that an attractor regime has been reached, where 
d T (j) = /(</>). In the context of Horndeski’s theory, it is known that on superhorizon scales the nonlinear curvature 
perturbation is conserved in this regime so long as the so-called decaying mode can be neglected, which implies that 
on superhorizon scales the nonlinear comoving curvature perturbation is both disformally invariant and conserved in 
any theory that is disformally related to Horndeski’s, provided d T cj> = /(</) and the decaying mode can be neglected. 
In addition, we found that the vector and tensor perturbation are invariant under a general disformal transformation 
at the nonlinear level. 

Finally, we discussed the healthiness of theories disformally related to Horndeski’s in relation to the presence of 
Ostrogradsky instabilities. Focusing on a toy model, we saw that the appearance of higher derivatives in an action 
that results from making a field transformation involving time derivatives of some of the fields is accompanied by 
the appearance of second class primary and secondary constraints in the Hamiltonian analysis. These constraints 
therefore remove the spurious additional degrees of freedom so that, as one would naively expect, the number of 
degrees of freedom is unaffected by the field redefinition. This is in agreement with the results of a recent analysis by 
Domenech et al. in [29], where it was shown that the number of degrees of freedom is indeed unaffected by disformal 
transformations, so long as the transformations are invertible. 

Note added: Our results are in agreement with those presented in Ref. [46], which appeared on arXiv at the same 
time as the present article and contains some material that overlaps with the nonlinear analysis we presented in § 111 . 
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Appendix A: Vanishing of e s in Horndeski’s theory 

In this appendix we derive the Poisson equation (2.29) in Horndeski’s theory, which is used in §11 to argue that 
under reasonable assumptions the comoving density perturbation e s , defined in ( 2 . 21 ), is suppressed on superhorizon 
scales. This result is a generalization of the result obtained in [27] for a subclass of Horndeski’s theory with G 4 = 1/2, 
G 5 = 0. 

Taking the action (2.24), assuming a flat FLRW background with line element of the form ds 2 = —N 2 {t)dt 2 + 
a 2 (t)Sijdx l dxi and varying the action with respect to N(t) gives us one of the background equations of motion [33] 

s = Y, £i = o, (Ai) 


where 


£2 = 2 XK x - K, (A 2 ) 

£.3 = 6 Xj)HG 3X - 2 XG 34 ,, (A3) 

£4 = -6 H 2 G a + 24 H 2 X{G 4X + XG 4X x) - 12 HX^G^x - 6 H^G 44> , (A4) 

£5 = 2H 3 X0{5G 5X + 2XG 5XX ) - 6ff 2 X(3G 50 + 2 XG^ x ). (A5) 


Whilst there are two additional background equations corresponding to a variation with respect to a(t) and we 
will not use them in the following analysis, so we omit them here. 

Turning next to linear perturbations, we take a line element of the form given in (2.8), but focus on the scalar 
perturbations. In addition, we fix the spatial gauge such that E = 0. The equations of motion for the perturbations 
are obtained by expanding the action (2.24) to second order and varying it with respect to each of the perturbation 
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variables. Here we will only need two of the equations - the two constraint equations obtained by varying the action 
with respect to A and B - and they take the form [33] 


k 2 ■ k 2 

—A 4 ip + A 4 A H—^ (—A 3 ip — aBA 5 ) = —A 2 8 (p — A 6 ^rS(p + /r \ 

or or 

-Cij> + C 3 A = -C 284 > - C 4 8 (p, 

where the coefficients are dependent on background quantities and are given as 

A 4 = 2E, 


60, 

2(E + 3i70) 

^2 = • » 

A 3 = 2Qt, 


<P 


- 20 , 

1 2(0 - HQ t ) 

A 6 = ■ , 

M = ^4>i 


<i> 

2 Qt, 

2(0 - HQ t ) 

^2 — • , 

£> 

11 

1 

to 

® 


0 


C 4 = -T- 

02 


2(H(j) - H<j))g T - 2 4>e 


with E, 0 and Qt being dehned as 


y = X X + l„X 0 = -i^ C =2 

X dx + 2 H dH 1 6 dH' Qt 


Ga — 2XGax — 


4 X — X x — G50) 


(A6) 

(A7) 


(A8) 


(A9) 


The easiest way to obtain (2.29) is to consider the unitary gauge, where Sep = 0. In this gauge we have ip = —TZ C , 
A = —e s /4> 2 and B = — (\& + lZ c )/(aH) and Eqs. (A 6 ) and (A7) thus reduce to 


k 2 


A 3 


a 4 tz c - a 4 ^ + — ( A 3 n c + -^(4- + 7 Z c )) = 0, 
C 4 ii c -C 3 ^- = 0 . 


(A10) 

(All) 


Eliminating 1Z C from these equations we arrive at (2.29). 

One can, of course, derive (2.29) without having to fix the gauge. To see this explicitly, we eliminate ip from 
Eqs. (A 6 ) and (A7) to obtain 

[a 2 - XXj Sep + U 4 - XXj A - ^ 6 <p=^ ( A 3 iP + aBA 5 - A 6 S<p) . (A12) 

Considering the left-hand side (l.h.s.) of (A12) first, it can be re-written as 


l.h.s. 


1 ( AlC 3 
cp 2 V Ci 



Using (A 8 ) one can then confirm that 


2 fA 2 C 1 -A 1 C 2 \ 

\a 4 c 3 -a 4 c 4 ) 


Sep - ep 2 A - </> 2 


AiC 4 + /rUi \ 

A&-A 4 CJ 


A 2 C\ — A 1 C 2 1 
A 1 C 3 — Ci A 4 ip 


(A13) 


(A14) 


In evaluating the coefficient of the 6 <p term we note that as a result of the background equation of motion £1 = 0 we 
have 


d£_ 

dt 


d£: 

dep^ 


d£ 


X 


d£ 


H = 0, 


(A15) 


dX dH 

from which we are able to obtain an expression for fi = Using this result, and noting that X = 2Xep/ep 1 we find 


Altogether we thus obtain 


AiC 4 fxC\ ep 

A1C3 — A 4 C'i (p 2 


l.h.s. 


1 

ep 2 


A 1 C 3 



e s - 


(A16) 


(A17) 
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Turning next to the right-hand side (r.h.s.) of (A12), we are able to re-write it as 


v fc2 

r.h.s. =-r 

a z 


A 3 n c + ^ (* + n c 


A 3 H , A s . „, 
—• b A 6 d f ) S(j) 


(A18) 


Using (A 8 ) one can then show that A 3 H + < j)A 3 + A 5 = 0. As such, we have recovered (2.29) without fixing the time 
slicing. 


Appendix B: Equivalence of equations of motion in disformally related frames 


In this appendix, following the analyses of [31, 40], we show that the gravitational equations given in (4.2) are 
equivalent to the gravitational equations of motion determined by varying (4.1) with respect to g so long as the 
disformal transformation relating g pl , and g^ is invertible, namely the invertibility condition ( 2 . 6 ) is satisfied. 

We will use the following relation between Tfif and , i.e. the inverse relation of (4.5), which is given as 

T% = a^a(a-2Xp) joT#' + ( ax g a p + 0 x d a d>ds<p) T*' 3 4 5 j . (Bl) 


Here, the coefficient in front of the large bracket on the right hand side does not vanish as (2.6) is satisfied. Using 
the above relation it is possible to rewrite (4.2) as 


« (V + + | (e% + \t£\ ( a x9pa + dxd p 0 d a <P) = 0. (B2) 

Contracting this equation once with g pu and once with d p 0d v (j) we obtain the two equations 

(a~Xa x )P-X(3 x Q = 0, (B3) 

2 X 2 a x P + (a + 2X 2 (3 x )Q = 0, (B4) 

where 

P = (stf + \n, V ) 9 ^ and Q = (stf + l -T^ d p cpd v 0 , (B5) 

The solution to these two equations is P = 0 and Q = 0, so long as the matrix M, given as 


. r _ fa — Xax —Xj3 x \ 

~ V 2A 2 « a - a + 2X 2 p x J 

is invertible. Substituting P = Q = 0 into (B2), and assuming a ^ 0, we thus recover 

cV-v 1 __ n 

°H +9 1 h — u ’ 


(B 6 ) 


(B7) 


i.e. we have recovered the gravitational equations obtained on minimising (4.1) with respect to g pu . The case when 
M is not invertible corresponds to 

a(a - Xa x + 2 X 2 /3 x ) = 0 , (B 8 ) 

meaning that we require a yb 0 and a — Xa x + 2X 2 (3 x 7 ^ 0 , which is the case so long as the invertibility condition 
( 2 . 6 ) is satisfied. 
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